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FIRST INTRINSIC ANISOTROPY OBSERVATIONS WITH THE COSMIC BACKGROUND IMAGER 
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ABSTRACT 

We present the first results of observations of the intrinsic anisotropy of the cosmic microwave background 
radiation with the Cosmic Background Imager from a site at 5080 m altitude in northern Chile. Our observations 
show a sharp decrease in C/ in the range / = 400-1500. The broadband amplitudes we have measured are <57b an d = 
58.7!^ fiK for / = 603^ and <5r band = 29.7^| M K for / = 1 190±ffJ, where these are half-power widths in /. Such 
a decrease in power at high I is one of the fundamental predictions of the standard cosmological model, and these 
are the first observations which cover a broad enough I range to show this decrease in a single experiment. The Q 
we have measured enables us to place limits on the density parameter, il tot < 0.4 or £l to t > 0.7 (90% confidence). 

Subject headings: cosmic microwave background — cosmology: observations 



1. INTRODUCTION 

In standard cosmologies, spatial temperature variations in the 
cosmic microwave background radiation (CMBR) are closely 
related to the primordial density fluctuations which gave rise 
to the formation of all structure in the universe (Peebles & Yu 
1970; Sunyaev & Zeldovich 1970). The angular power spec- 
trum of these temperature variations on the celestial sphere, C;, 
yields a direct estimate of the prime cosmological parameters 
and provides a fundamental link between particle physics and 
cosmology (e.g., Kamionkowski & Kosowsky 1999). Since ra- 
dio interferometers sample the angular power spectrum directly 
and are straightforward to calibrate, they provide a simple and 
direct determination of C/. Here we report the first observations 
of the CMBR with the Cosmic Background Imager (CBI). 

2. THE COSMIC BACKGROUND IMAGER 

The CBI is a radio interferometer with 13 0.9 m diameter an- 
tennas mounted on a 6 m tracking platform. It operates in 10 
1 GHz frequency channels from 26 to 36 GHz. The instanta- 
neous field of view and the maximum resolution are ~ 45' and 
~ 3' (FWHM). The instrument has an altitude-azimuth mount, 
and the antenna platform can also be rotated about the optical 
axis to increase the aperture-plane (u,v) coverage and to facil- 
itate polarization observations. The antennas have low-noise 
broadband high-electron mobility transistor (HEMT) amplifier 
receivers with ~ 25 K noise temperatures. The typical system 
noise temperature averaged over all 10 bands is ~ 30 K, in- 

1 See http://www.astro.caltech.edu/~tjp/CBI. 

2 See the MAP Web site at http://map.gsfc.nasa.gov/. 



eluding ground spillover and atmosphere. The frequency of 
operation of the CBI was chosen as a compromise between 
the effects of astronomical foregrounds, atmospheric emission, 
and the sensitivity that can be achieved with HEMT amplifiers. 
Details of the instrument design may be found in Padin et al. 
(2000a,b, S. Padin et al. 2000, in preparation) and on the CBI 
web site. 1 The CBI is located at an altitude of 5080 m near 
Cerro Chajnantor in the Atacama desert in northern Chile. This 
site was chosen because the atmospheric opacity is low and the 
CBI can operate at the thermal noise limit much of the time. 
The instrument was assembled and tested on the Caltech cam- 
pus during 1998 and 1999 and shipped to Chile in 1999 August. 
Installation of the telescope and site infrastructure were com- 
pleted by the end of 1999, and the full instrument has been in 
operation since early 2000 January. 

The CBI is sensitive to multipoles in the range 400 < / < 
4250, where these values reflect the half-power widths of the 
window functions on the shortest and longest baselines. The 
CBI complements BOOMERANG (de Bernardis et al. 2000; 
Lange et al. 2001), the Degree Angular Scale Inteferometer 
(DASI) (Halverson et al. 1998), Microwave Anisotropy Probe, 2 
MAXIMA (Hanany et al. 2000; Balbi et al. 2000), and the 
Very Small Array (Jones & Scott 1998), which cover the range 
50 < / < 1000. DASI is a sister project to the CBI, and the CBI 
and DASI designs were chosen to complement each other. The 
CBI control software and correlator and receiver control elec- 
tronics were duplicated by the DASI team for the DASI project. 
Together these two interferometers cover the multipole range 
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100 < Z < 4250. 

3. OBSERVATIONS 

The antenna platform of the CBI permits a wide variety of 
antenna configurations. For observations during the test phase 
(2000 January-April) we chose a configuration with the anten- 
nas around the perimeter of the platform, which provided easy 
access to the receivers and fairly uniform («,v) coverage, en- 
abling us to test the full range of CBI baselines. We report 
here only observations on baselines corresponding to / < 1510, 
which account for 25% of the data in this ring configuration. 
We do not report on the higher multipole bins because these are 
more dependent on the bright-source and statistical faint-source 
corrections, which are still preliminary. These results, and mo- 
saicked observations which significantly increase the resolution 
in I, will be presented elsewhere. 

We based our flux density scale on observations of Jupiter, 
assuming 7j up i tel - = 152 K at 32 GHz, with 5% uncertainty (Ma- 
son et al. 1999). The spectral index of Jupiter is not constant 
between 26 and 36 GHz (Wrixon et al. 1971), so we used Tau- 
rus A as our prime calibrator. Taurus A is slightly resolved 
with the CBI, but it can be well fitted by an elliptical Gaussian 
model. We referenced the 32 GHz flux density of Taurus A to 
that of Jupiter and transferred this to the other frequency chan- 
nels assuming a-r. dums A = -0.30, where S oc v a (Mezger et al. 
1986). Any uncertainties resulting from this extrapolation are 
< 1%. 

For our intrinsic anisotropy observations, we selected regions 
at Galactic latitudes above 20°, for which the synchrotron emis- 
sion and the IRAS 100 [im emission are relatively low, and 
which avoided bright point sources detected in the NRAO Very 
Large Array 1.4 GHz Sky Survey (NVSS) (Condon et al. 1998). 
These regions are also in the declination range -5° < 5 < -2°, 
which permits point-source monitoring with the Owens Valley 
Radio Observatory (OVRO) 40 m telescope (see § 5). 

Daytime observations of the CMBR are not possible because 
radio emission from the Sun contaminates the visibilities. Ra- 
dio emission from the Moon causes similar problems, and in- 
trinsic anisotropy observations are restricted to fields more than 
60° from the moon. The weather during 2000 January-April 
was uncommonly poor, and we lost 50% of the nights as a 
result of bad weather. In the remaining nights the observing 
conditions were superb, and the sensitivity was limited by the 
system noise. 

The dominant systematic contamination in the CBI observa- 
tions is due to ground spillover. The level of contamination on 
the 1 m baselines is typically between a few tens and a few 
hundreds of mJy; in rare instances it can be as high as a few 
Jy, and it falls off with increasing baseline length. The ground 
signals can be recognized in the maps because they are not con- 
fined to the primary beam area, as celestial signals are. We 
have successfully removed the ground signal by differencing 
visibilities measured on two fields at the same declination sep- 
arated by 08™ in right ascension. The observations were alter- 
nated between two fields on an 8 minute timescale, so that both 
fields were observed at the same position relative to the ground. 
In addition to filtering out ground spillover, this differencing 
strongly rejects crosstalk and other spurious instrumental sig- 
nals. As discussed below, any residual spurious signals amount 

1 /2 

to < 1.3% of C,' and may therefore safely be ignored. The 
differenced observations reported here are between two fields 
centered at 08 h 44 m 40 s , -03° 10' and 08 h 52 m 40 s , -03° 10'; and 



two fields centered at 14 h 42 m , -03° 50' and 14 h 50 m , -03° 50' 
(J2000). To search for residual spurious signals in excess of the 
thermal noise level in our differenced visibilities, we needed an 
accurate estimate of the thermal noise level. We obtained this 
estimate by computing the rms of successive 8.4 s integrations 
in the differenced visibilities. The typical noise level computed 
in this way was 2 Jy, in agreement with the (less accurate) noise 
level determined from the system temperature. The means of 
the 8™ scans were much smaller than 2 Jy, so this method pro- 
vided an accurate estimate of the thermal noise level. To test for 
residual spurious signals, we divided our differenced visibilities 
into two parts in three ways: (1) pre- versus posttransit observa- 
tions; (2) middle half versus outer quarters of hour-angle range; 
and (3) first versus second half of epochs on each field (pre- 
March 23 vs. post-March 23). 

In each case, we subtracted the singly differenced visibilities 
point by point in the («,v) plane. In cases 1 and 2 the doubly 
differenced visibilities were consistent with the expected noise. 
In case 3 there was an excess amounting to a 1.3% contamina- 
tion in C 1 / 2 . 

4. ANALYSIS AND RESULTS 

We observed for 58.5 hr on each of the 08 h fields and for 
16.15 hr on each of the 14 h fields. The sky signal is clearly vis- 
ible in the differenced images on the 100 and 104 cm baselines 
(Fig. 1). 

The extraction of the angular spectrum from visibility mea- 
surements is straightforward (White et al. 1999). The covari- 
ance matrix of the observations is the matrix of the covariances 
between all the visibility measurements: 

C = M + N, (1) 

where M and N are the sky and noise covariance matrices. We 
assume that the noise on different baselines and at different fre- 
quencies is uncorrelated, i.e., that N is diagonal. The sky co- 
variance matrix is 

M jk =(V(vL j ,u j )V*(n k ,iy k )) (2) 




d 2 \A(uj - v, Vj)A*(u k - v, v k )S(v, v b v k ) (3) 



for two visibility points j, k, where A(u, v) is the Fourier trans- 
form of the primary beam at frequency v, u is the baseline vec- 
tor in wavelengths, and S(\, Vj, v k ) is a generalized power spec- 
trum of the intensity fluctuations (Hobson et al. 1995). The 
effective weighting of this power spectrum defines the window 
function. The indices, j and k, run from 1 to n where n is the 
number of distinct (u, v) points. We average all the data taken at 
different times for each (m, v) point before doing the maximum 
likelihood calculation. 

The generalized power spectrum is related to C/ by 

S(v,Uj,u k )=(^^j v)vlz(v^g(v k )C u (4) 

where l+\ = 2tt\\\ and the g factor is a small correction for the 
difference between the Rayleigh-Jeans and Planck functions. 
We can test a hypothetical power spectrum, [Q], by forming 
the likelihood function 

UlCii) = — ex V {-V*(Uj)C] k l V(Uk)}. (5) 

The cross-correlation between the signals received from two 
fields separated by 08 m in R.A. is negligible, so the expected 
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Residual LL map. Array: CBI 
C0844-0310 at 31.000 GHz 2000 Jan 12 



Dirty LL beam. Array: CBI 
C0844-0310 at 31.000 GHz 2000 Jan 12 



Clean LL map. Array: CBI 
C0844-0310 at 31.000 GHz 2000 Jon 12 





Right Ascension (orcmin) 

Map center; RA; 08 44 40.000. Dec: -03 09 59.995 (2000.0) 
Displayed range: -0.056 -» 0.0606 Jy/beam 



Displayed range: -0.386 ■* 1 



Right Ascension (arcmin) 
Map center: RA: 08 44 40.000, Dec: -03 09 59.995 (2000.0) 
Map peak: -0.033 Jy/beam 
Beam FWHM: 16.4 x 16.3 (arcmin) at -49.2° 



FIG. 1. — Differenced image of the 08 h field observed on the nine 100 cm and 104 cm CBI baselines in the test configuration (left), the corresponding point-spread 
function (center), and an image deconvolved with the Hogbom CLEAN algorithm (right), assuming that the emission is confined to the primary-beam area. The 
10 frequency channels have been combined, and the synthesized beamwidth ~ 16' (FWHM). The circle at radius 44' indicates the primary-beam area which was 
cleaned. The structures outside the circle in the uncleaned image (left) are due to sidelobes. These structures disappear after cleaning within the primary-beam area, 
indicating that they are caused by sidelobes. 



variance of the differenced visibilities is twice the variance of 
the undifferenced visibilities. Our parametric model for C/ con- 
sists of two parameters, these being the amplitudes of Q in the 
two ranges / < 900 and / > 900, assuming 1(1 + l)Q is con- 
stant in each range. In the test configuration there is a gap 
in our (u,v) coverage between / = 800 and / = 1000 which 
makes this a natural division. The band-power window func- 
tions, here approximated by sums of the single-baseline single- 
channel window functions within each / bin, may be charac- 
terized by I = 603!||> and / = 1190^ (half-power widths). 
The maximum likelihood broadband signals we measure are 
<57band = [1(1+ 1)C / /(2tt)] 1 /2 x T cmb = 58.7!" /zK for the / = 603 
bin and (57b an d = 29.7^ | for the I = 1 190 bin. The error bars 
indicate the points at which the likelihood has dropped by e" 5 
(which are within 10% of the 68% integrated probability val- 
ues). The results for the two fields, shown in Figure 2, agree to 
within the uncertainties. For the lower I bin, the uncertainties 
are dominated by sample variance, and they would be decreased 
by < 2% in the absence of thermal noise. For the upper / bin, 
the uncertainties would be decreased by 3 1 % and 54% for the 
08 h and 14 h fields, respectively, in the absence of thermal noise. 

Our maximum likelihood analysis has been tested using soft- 
ware written independently by two of the authors, with no com- 
mon code between the packages and using significantly differ- 
ent implementations for important steps, such as evaluation of 
the window function, binning and maximization algorithm. We 
have analyzed both the real data and simulated differenced data 
sets generated by realizations of known power spectra together 
with realistic noise and point sources. The results from these 
two software packages are in excellent agreement, and the sim- 
ulations recover the original input power spectra. We have also 
generated 54 simulations of differenced sky images based on 
our observed band powers in the two bins (Fig. 2) and com- 
pared the rms signal, measured within the primary-beam area 
in these simulations, with the rms fluctuations in the primary 
beam measured in actual observations of 54 differenced fields, 
to be published elsewhere. Both the means and the distribu- 



tions of the rms values for the observed and simulated fields are 
in excellent agreement. Thus, we are confident that our derived 
spectrum is a reliable representation of the signal that we have 
detected on the sky. 

5. FOREGROUNDS 

Radio galaxies and radio-loud quasars are a source of con- 
fusion at CBI frequencies and angular scales, so we equipped 
the OVRO 40 m telescope with a four-channel 26-34 GHz 
receiver for point-source monitoring, and observed all of the 
sources in the NVSS with Si.4gh z > 6 mJy in our CBI fields. 
Those sources detected with the 40 m telescope at the 3 a level 
(S30GHZ > 6 mJy) have been subtracted from our CBI visibility 
data using the flux densities measured on the 40 m telescope. 
This reduced the levels of <57b an d measured in the lower and up- 
per I bins by 0.5% and 1%. We have also applied corrections 
based on the source count statistics (White et al. 1997) to ac- 
count for point sources with S30GHZ < 6 mJy, which have not 
been subtracted individually from our visibility data. The cor- 
rections amounted to decreases in <57b an d of 1.6% and 8.3% in 
the lower and upper I bins. These corrections have been applied 
to the band powers given in § 4. The uncertainty in the statisti- 
cal correction, ~ 20%, makes <C 0.1 K difference to the errors 
in both bins. 

It is unlikely that diffuse Galactic foreground emission is a 
significant contaminant in our observations. The expected rms 
fluctuations due to Galactic synchrotron emission on angular 
scales 5-30' are < 9 fiK (e.g., Tegmark et al. 2000). In the 
RING5M experiment, Leitch et al. (1997) detected Galactic 
emission at 14.5 and 32 GHz with a spectrum consistent with 
free-free radiation but a much higher level than predicted from 
Ha measurements. We have therefore made 14.5 GHz obser- 
vations with the OVRO 40 m telescope along a strip at decli- 
nation -5° over the right ascension range h -24 h . The beam 
and beamthrow, 7.'4 and 22'.2, are fairly well matched to the 
CBI angular scales in the lower I bin, so, after correcting for 
the window function, we may use these observations to esti- 
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FIG. 2. — CMBR anisotropy spectrum determined from CBI observations. The triangles and squares show results on the 08 h and 14 h differenced fields; the circles 
show the results of a joint maximum likelihood analysis of both differenced fields. The individual 08 h and 14 h field results are offset in / for clarity. The window 
functions for each bin are shown as dashed lines. The solid curve represents a flat model universe with Hq = 75 km s~' Mpc~' , Q\,h 2 = 0.019, and Q ct i m = 0.2. 



mate the possible level of contamination in our CBI observa- 
tions. If all of the signal seen at 14.5 GHz at \b n \ > 5° is at- 
tributed to anomalous emission with the same spectral proper- 
ties as seen in the RING5m data, then this amounts to a compo- 
nent in our CBI observations (57t, an d = 20 /iK due to anomalous 
foregrounds. Subtraction in quadrature from the signal we have 
detected in the first I bin would reduce our observed <57b an d by 
7%. 

We have measured the temperature spectral index, (3 = 
\n{T\/T2)/\Yi{vi/v2), with 1 month of data from the CBI in a 
more compact configuration, optimized to measure both the an- 
gular spectrum and the radio-frequency spectrum of the CMBR. 
We find that the signal is primarily CMBR and not Galactic. We 
used a maximum likelihood analysis with f3 as a free parame- 
ter to determine that /? = 0.0±0.4 (la error) in the lower I 
bin. If as much as 21% of the <57b an d in this I bin were due to 
a free-free foreground component with spectral index (i = -2.1, 
while the remainder was CMBR with f3 = 0, the spectral index 
measured would be < -0.8, which is ruled out at the 2 a level. 
A 15% synchrotron foreground component with spectral index 
(i = -2.7 can be ruled out at the same level. 

6. DISCUSSION 

A decrease in C/ at high /, caused by photon diffusion and 
the thickness of the last scattering region, is a fundamental 
prediction of the standard cosmological model (Silk 1968), 
and this is the first time that such a decrease has been de- 
tected in a single experiment. It is also the first time that 
anisotropy has been detected at / > 1000. The levels of <57b an d 
detected with the CBI are consistent with observations at high 
/ made over the last 12 years (Readhead et al. 1989; Scott et 
al. 1996; Church et al. 1997; Baker et al. 1999; Leitch et al. 
2000; Holzapfel et al. 2000; Subrahmanyan et al. 2000). The 
level of <57t, an d we measure at / ~ 600 is a factor 1.5 higher 
than that found by BOOMERANG and a factor 1.3 higher 
than that found by MAXIMA; here we have used the best fit 
spectrum to the BOOMERANG+MAXIMA+DMR data (Jaffe 
et al. 2001) to extrapolate the BOOMERANG and MAXIMA 
data. The additional power <57j^ nd detected by the CBI is sig- 
nificant at the 1.8 a level (BOOMERANG) and the 1.4 a level 
(MAXIMA), where a includes calibration uncertainties of 10% 
(BOOMERANG), 4% (MAXIMA) and 5% (CBI), pointing un- 
certainties of 11% (BOOMERANG) and 5% (MAXIMA), the 



uncertainty of 13% in the (57t, an d measured in the lower I bin 
(CBI), and an estimated uncertainty of 8% in the <57b an d mea- 
sured by BOOMERANG and MAXIMA in the multipole range 
300 < I < 700. It is important to determine whether these dif- 
ferences between the CBI and BOOMERANG-MAXIMA are 
real. The RING5M experiment (Leitch et al. 2000) reported 
<57band = 59±^K at / ~ 600, which agrees well with the CBI 
value and is discrepant at the 1.8cr level with BOOMERANG 
and at the 1 A a level with MAXIMA. The CAT values (Scott et 
al. 1996; Baker et al. 1999) are intermediate between the CBI 
and BOOMERANG and MAXIMA values, but any differences 
are significant only at the ~ 1 a level. 

We have used the likelihoods of our data to explore limits 
on the cosmological parameters shown in Table 1 for both flat 
and open model universes with power-law density fluctuation 
spectra having slope n = 1 using CMBFAST (Seljak & Zaldar- 
riaga 1996). £1 is the density parameter, and subscripts "tot," 
"b," "m," "cdm," and "A" refer to the total, baryonic, matter, 
cold dark matter, and cosmological constant contributions to 
the density parameter. The ranges and intervals of the param- 
eters are shown in Table 1. For the flat models we find the 
likelihood peaks at £l\,h 2 = 0.009 and drops by a factor of 2 
at n b h 2 = 0.019 (e.g., Buries & Tytler 1998; O'Meara et al. 
2001), and by a factor of 3 at il\,h 2 = 0.03. For the open models 
we have assumed uniform priors for Ho, flbh 2 , Q c dm, and f2 to t, 
with £l m being uniformly distributed between 0.2 and £! tot , over 
the ranges indicated in Table 1, and we find that fl tol < 0.4 or 
f2 to t > 0.7 at the 90% confidence level. 
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Table 1 

Parameter Space for Cosmological Likelihood Analysis 



Flat Models Open Models 



n b h z 


0.003^0.03(0.003) 


0.01 


-> 0.03 (0.01) 


^cdm 


0.1 0.5(0.04) 






£l m 




0.2- 


-fit* (0.1) 


ttA 








^tot 


1 


0.2- 


-> 1.0(0.1) 


h 


0.50^0.90(0.04) 


0.50 


^0.80(0.05) 



The step size for independent variables is shown 
in parentheses. 
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